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ARID LAND GEOGRAPHY 


基于 反射 率 和 党 度 温 度 的 锡林郭勒 积 雪 深度 估算 


(锡林浩特 国家 气候 观象台 ,内 蒙古 锡林浩特 026000) 


摘 要 : 积 雪 是 我 国 西北 


干旱 半 干 旱 区 重要 的 水 资源 ,也 是 影响 全 球 气 候 变 化 的 重要 因子 之 一 。 


目前 光学 影像 反射 率 和 雷达 亮 温 数据 是 积 雪 逮 感 领 域 的 主要 数据 ,本 文 首 次 结合 两 类 遥感 数据 估 


算 积 雪 深 度 ,并 比较 偏 最 小 二 乘法 和 机 器 学 习 算 法 (人 工 神 经 网 络 、 支 持 向 量 机 和 随机 森林 算法 ) 
在 积 雪 深度 估算 方面 的 表现 。 以 锡林郭勒 嚼 2012 一 2015 年 积 雪 深度 数据 为 例 ,基于 反射 率 和 亮度 


温度 相 结合 的 积 雪 深度 估算 精度 优 于 单个 数据 源 , 且 随机 森林 算法 表现 最 好 , 均 方 根 误差 为 2.93 
cm, 满足 实 际 应 用 的 需求 。 研 究 结果 对 我 国 西北 地 区 水 资源 分 布 . 生 态 环境 评估 等 研究 具有 重要 


意义 。 
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职 雪 是 影响 全 球 气 候 变 化 的 重要 因子 之 一 , 它 
对 地 气 系 统 辐射 平衡 的 影响 早已 被 世界 各 国 科学 
家 证 实 。 积 雪 也 是 我 国 西北 半 王 旱地 区 极为 重要 
的 水 资源 ,在 区 域 水 分 平衡 中 发 挥 着 极为 重要 的 作 
用 ,与 农 牧 业 生 产 有 着 十 分 密切 的 关系 "。 冬 季 积 
雪 丰 厚 可 以 对 农作物 和 牧草 起 到 保温 保 丧 的 作用 ， 
为 春季 农业 生产 提供 良好 的 水 资源 保障 ”。 因 此 ， 
积 雪 一 直 是 人 们 气候 变化 研究 的 重点 内 容 之 一 , 安 
W .科学 的 评估 积 雪 资源 状况 对 于 水 资源 循环 和 全 
球 能 量 交 换 人 研究 具有 重要 意义 。 积 雪 深 度 作 为 
积 雪 资 源 评估 中 的 重要 参数 之 一 ,如 何 准确 地 估算 
职 雪 深度 ,一直 是 积 雪 遥感 研究 中 的 重要 内 容 5= 。 

目前 ,估算 积 雪 深 度 的 遥感 数据 可 以 分 为 光学 
遥感 数据 和 被 动 微波 遥感 数据 。 在 光学 遥感 领域 ， 
陈 乾 “和 周 咏 梅 "分 别 建立 了 祁连山 区 和 青 南 高 原 
积 雪 深度 与 AVHRR 反射 率 的 统计 关系 。 近 年 来 ， 
具有 更 高 空间 分 辨 率 和 光谱 分 辩 率 的 MODIS 数据 
thie BAF REET, TER oh 
H, CHANG 等 "在 辐射 传输 和 米 氏 散射 的 理论 基础 
上 建立 了 积 雪 深度 与 19 GHz 和 37 GHz 水 平 极 化 的 
亮 温 差 的 统计 方法 。 于 惠 等 、 卢 新 玉 等 "中 也 
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采用 了 类 似 的 研究 思路 估算 了 中 国 地 区 的 积 雪 深 
E, PAWE MEDEE RRRA T RE 
FY-3B 和 FY-3C MWRI 微 波 辐射 数据 在 雪 深 估算 的 
潜力 ,取得 了 较 好 的 估算 效果 。 经 过 不 断 的 发 展 与 
改进 ,基于 亮 温 的 方法 逐渐 成 为 了 被 动 微波 遥感 估 
算 积 雪 深 度 应 用 最 广泛 的 算法 之 一 。 

在 积 雪 深度 估算 领域 ,光学 遥感 无 法 穿 透 积 雪 
表层 获得 积 雪 内 部 信息 , 当 积 雪 深 度 超过 一 定 范 
FB] , 某 些 波段 的 反射 率 达 到 饱和 ,限制 了 估算 精 
度 。 与 光学 遥感 相 比 ,被 动 微 波 遥 感 在 穿 透 能 力 与 
全 天 候 工作 特性 方面 具有 先天 性 优势 ,可 用 于 进行 
全 天 候 的 积 雪 监 测 王 ,更 加 适用 于 积 雪 深度 估算 ， 
但 是 其 空间 分 辨 率 较 低 ,通常 在 25 km 至 100 km, Më 
以 满足 局 地 尺度 空间 分 辩 率 研究 的 需求 。 已 有 研 
究 综合 利用 光学 遥感 数据 和 被 动 微波 遥感 数据 进 
行 积 雪 识别 “” ,然而 综合 利用 两 类 数据 监测 积 雪 
深度 的 研究 还 未 见报 道 。 因 此 ,本 文 的 首要 研究 目 
的 就 是 评估 光学 遥感 和 被 动 微波 遥感 数据 结合 在 
积 雪 深度 估算 的 潜力 。 

考虑 到 积 雪 深度 估算 的 复杂 性 ,简单 的 线性 模 
型 有 可 能 使 得 估算 结 果 存 在 很 大 的 不 确定 性 "|。 
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偏 最 小 二 乘法 ”可 以 克服 自 变 量 间 的 共 线 性 ,在 


率 数据 。 该 数据 最 大 限度 的 去 除了 云 的 影响 ,保留 


遥感 领域 具有 重要 应 用 。 除 此 之 外 ,人 工 神经 网 
络 ”支持 向 量 机 和 随机 森林 算法 “等 机 带 学 
习 算 法 也 吸引 了 众多 遥感 学 者 的 关注 ,其 中 ,神经 
网 络 可 以 处 理 复杂 的 映射 关系 ,支持 向 量 机 可 以 将 
低 维 数据 转 到 高 维 空间 进行 回归 ,随机 森林 算法 学 
习 过 程 快 .和 鲁 棒 性 强 。 因 此 ,本文 的 第 二 个 研究 目 
的 就 是 比较 4 种 方法 积 雪 深度 的 估算 精度 ,为 北方 
地 区 积 雪 深度 遥感 估算 提供 理论 参考 依据 。 


1 研究 区 概况 


锡林郭勒 盟 (112° ~ 120° E,42.5° ~ 46.7°N) 位 
于 内 蒙古 自治 区 中 部 ,海拔 900 ~ 1 300 m, 属 于 温带 
干旱 半 干 旱 大 陆 性 气候 ,冬季 寒冷 ,夏季 较 热 ,年 均 
气温 1~2%。 全 盟 大 部 分 地 区 积 雪 时 期 集中 在 11 
月 ~ 次 年 3 月 份 。 


2 数据 与 方法 


2.1 积 雪 深度 数据 

本 文 积 雪 深度 数据 来 源 于 锡 盟 15 个 气象 站 点 
的 实地 观测 数据 ,时间 覆 盖 范 围 为 2012 一 2015 年 ， 
观测 的 积 雪 深度 最 大 值 为 37 cm。 根 据 MODIS 地 
物 分 类 图 可 知 , 这 些 气 象 站 点 的 下 垫 面 均 为 草地 
(图 1)。 
2.2 MODIS 反射 率 数据 

本 文 下 载 了 『 MODIS/Terra 提供 的 8 d 合成 反射 
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1 研究 区 地 理 位 置 及 气象 站 点 分 布 
Fig. 1 Location of the study area and distribution of the 


meteorological stations 


了 质量 较 高 反射 率 数据 。 归 一 化 积 雪 指数 (Normal- 
ized Difference Snow Index, NDSI ) 利 用 积 雪 在 可 见 
光 高 反射 和 短波 红外 较 强 吸收 的 特点 ,和 常 被 用 于 区 
分 积 雪 与 非 积 雪 目 标 , 是 积 雪 遥感 领域 应 用 最 广 的 

BRS eT ,计算 公式 如 下 : 

NDSI=(R4 — R5)/(R4 - R5) (1) 

式 中 :R4 和 R6 分 别 为 MODIS 第 4 波段 和 第 6 波段 的 
反射 率 。 

2.3 FY-3B 亮 温 数 据 
风云 3B(FY-3B) 卫 星 是 我 国 发 射 的 第 2 代 极 轨 

气象 卫星 ,其 搭载 的 微波 成 像 仪 (MWRI) 提 供 了 5 个 
中 心 频 率 (10.65 .18.7.23.8 36.5 89.0 GHz) .2 种 极 
化 方式 (H 和 V) 的 亮 温 数据 ,总 共 10 个 波段 
MWRI 数 据 包 含 升 轨 和 降 轨 数据 ,过 境 时 间 分 别 为 
14:00 和 02:00 左 右 , 其 中 ,白天 升 轨 数据 用 于 本 文 
的 研究 。 

2.4 数据 处 理 流 程 

首先 ,根据 气象 观测 站 点 的 经 纬度 提取 与 积 雪 
深度 观测 时 间 对 应 的 MODIS 反 射 率 数 据 和 FY-3B 
亮 温 数据 ;其 次 ,计算 NDS7, 保 留 NDST 大 于 0.4 的 数 
据 ;最 后 ,反射 率 和 亮 温 数据 作为 回归 模型 的 因 变 
量 , 积 雪 深 度 作为 回归 模型 的 自 变 量 。 在 精度 验证 
阶段 ,本 文采 用 十 字 交 叉 验 证 方法 对 模型 进行 
评价 。 


3 ”结果 与 分 析 


3.1 积 雪 数据 概况 

基于 15 个 气象 观测 站 点 长 时 间 雪 深 数据 ,本 文 
总 共 使 用 了 256 个 样本 数据 ,其 中 积 雪 深 度 的 频率 
直方 图 见 图 2, 积 雪 深 度 观测 值 并 不 服从 正 态 分 布 ， 
不 满足 简单 线性 回归 正 态 分 布 的 要 求 , 这 也 说 明了 
简单 线性 回归 模型 不 适用 于 积 雪 深 度 的 反 演 。 随 
着 积 雪 深度 的 增 大 ,频率 降低 ,大 部 分 积 雪 深 度 观 
测 值 落 在 1 ~ 20 cm 范围 内 , 占 总 样本 个 数 的 86%。 
3.2 反 演 精度 比较 

表 1 汇 总 了 4 种 回归 方法 利用 光学 遥感 数据 、 
被 动 微波 遥感 数据 以 及 两 者 合成 数据 的 积 雪 深 度 
估算 精 度 。 偏 最 小 二 乘法 需要 确定 主 成 分 的 个 数 ， 
随机 森林 算法 需要 事先 确定 决策 树 的 变量 个 数 ,为 
了 避免 变量 个 数 对 偏 最 小 二 乘法 和 随机 森林 算法 
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1 5 9 13 17 21 25 29 33 37 
积 雪 深度 / cm 


图 2 积 雪 深 度 频率 统计 直方 图 


Fig.2 Frequency statistics histogram of snow depth 


表 1 四 类 回归 算法 在 不 同 遥感 数据 源 的 积 雪 
深度 反 演 精度 (RMSE) 


Tab.1 Inversion accuracy of four regression algorithm 


based on different remote sensing data /cm 
司 妥 算法 下 _ 
反射 率 亮 温 反射 率 + 亮 温 
偏 最 小 二 乘 8.63 8.60 7.12 
人 工 神经 网 络 6.88 5.53 3.92 
文 持 向 量 机 7.81 4.07 3.31 
随机 森林 算法 3.65 3.12 2.93 


注 :RMSE, 均 方 根 误差 


精度 影响 ,因此 表 1 中 的 精度 为 所 有 可 能 变量 情况 
下 的 平均 值 , 这 样 能 够 更 好 地 比较 不 同 回归 方法 以 
及 数据 源 之 间 的 精度 差异 。 

仅 使 用 反射 率 数 据 情况 下 ,随机 森林 算法 取得 
最 好 的 精度 (RMSE=3.65 cm) ,其 次 为 神经 网 络 和 支 
持 向 量 机 。 仅 使 用 亮 温 数据 情况 下 ,随机 森林 算法 
仍然 表现 最 好 (RMSE=3.12 cm) ,但 与 人 工 神 经 网 络 
相差 不 明显 。 偏 最 小 二 乘法 的 表现 始终 最 差 。 值 
得 注意 的 是 对 于 4 种 方法 ,使 用 亮 温 数据 佑 算 积 雪 
深度 的 精度 均 高 于 使 用 反射 率 数 据 , 这 可 能 归 因 于 


估算 雪 深 /cm 


= 
© 


ole 0 
0 10 20 30 40 0 10 20 30 
实测 雪 深 /cm 实测 雪 深 /cm 
(a) 偏 最 小 二 乘法 (b) 人 工 神经 网 络 


RMSE-392cm 


反射 率 数据 易于 饱和 。 

综合 反射 率 和 亮 温 数据 的 情况 下 ,表现 最 好 的 
仍然 是 随机 森林 算法 ,但 是 3 个 机 器 学 习 算 法 的 精 
度 差异 不 大 ,RMSE 为 2.93 ~ 3.92 em, 远 远 低 于 偏 最 
小 二 乘 方法 (RMSE=7.12 cm)。 偏 最 小 二 乘法 虽然 
有 很 广泛 的 应 用 ,但 结果 表明 ,该 方法 不 适合 于 积 
雪 深 度 的 估算 研究 ,而 机 器 学 习 算 法 能 够 比较 精确 
地 估算 积 雪 深度 。 

图 3 为 4 种 回归 方法 的 1:1 图 ,从 图 3a 中 可 以 明 
显 看 出 偏 最 小 二 乘法 估算 的 积 雪 深度 均 小 于 20 
cm, 佑 算 精 度 较 差 。 人 工 神经 网 络 的 实测 值 和 估算 
值 总 体 上 落 在 1:1 线 附近 ,但 是 较为 分 散 ( 图 3b)。 
对 支持 向 量 机 来 说 ,大 部 分 数据 点 落 在 1:1 线 上 ,但 
是 由 于 积 雪 深度 佑 算 的 复杂 性 ,有 一 些 点 偏离 1:1 
线 较 远 ,降低 了 支持 向 量 机 的 精度 (图 3c)。 表 现 最 
好 的 回归 算法 是 随机 森林 算法 (图 3d) ,大 体 上 所 有 
数据 点 均 较 为 紧凑 的 落 在 1:1 附 近 。 值 得 注意 的 
是 ,对 于 4 种 方法 来 说 ,估算 结果 在 积 雪 深 度 较 低 的 
时 候 出 现 高 估 ,而 在 高 积 雪 深 度 的 时 候 出 现 低估 。 

另外 ,对 于 随机 森林 算法 来 说 , 仅 使 用 亮 温 数 
据 与 综合 数据 的 精度 相差 0.19 cm ,也 就 是 说 单独 使 
用 亮 温 数据 已 经 可 以 取得 较 好 的 估算 精度 。 通 常 
情况 下 ,被 动 微波 遥感 数据 的 空间 分 辩 率 较 低 , 精 
度 高 ,而 光学 遥感 数据 的 空间 分 辩 率 高 。 融 合 两 类 
遥感 数据 有 望 获得 高 空间 分 辨 率 和 高 精度 的 积 雪 
深度 空间 分 布 图 。 

33 积 雪 深 度 空间 分 布 图 

本 文选 取 了 2016 年 1 月 4 号 锡林郭勒 盟 地 区 天 
气 晴朗 、 云 量 较 少 的 MODIS 影 像 数 据 ,图 4 为 MO- 
DIS 数据 500 mm 分 辩 率 的 真 彩色 图 , 锡 盟 绝 大 部 分 地 
区 被 积 雪 覆盖 。 从 VDS7 数 据 中 也 可 以 看 出 ,大 部 分 
地 区 NDSI 接 近 于 1( 图 5)。 


RMSE=2.93 cm A 


0 10 20 30 40 0 10 20 30 40 
实测 雪 深 /em 实测 雪 深 /em 
(c) 支持 向 量 机 (d) 随机 森林 算法 


图 3 实测 与 预测 积 雪 深度 的 1:1 图 
Fig. 3 The 1:1 plot of measured and predicted snow depth 
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图 4 锡林郭勒 盟 地 区 的 真 彩色 合成 图 
Fig.4 True color image of Xilin Col League 


图 6 为 融合 反射 率 和 亮度 温度 估算 的 积 雪 深度 
空间 分 布 图 (以 随机 森林 算法 为 例 , 其 他 算法 的 结 
果 未 显示 )。 锡 林 郭 勒 盟 地 区 东北 部 地 区 ( 东 乌 珠 
穆 沁 旗 和 乌拉 盖 ) 的 降雪 量 高 于 西部 地 区 (二 连 浩 
特 、 eal Slee allen 与 图 6 的 积 雪 深 度 空间 分 布 
图 一 致 。 最 高 积 雪 深 度 出 现在 东 乌 珠 穆 泥 旗 地 区 ， 
[se ne cmo 


4 讨论 与 结论 


反射 率 数 据 提 供 了 积 雪 


近 表面 的 光学 信息 ,而 


低 : -0.76 高 :1 


图 5 锡林郭勒 盟 地 区 MODIS-MDS7 
Fig. 5 MODIS-NDST image of Xilin Gol League 


积 雪 深度 /em 


D“ 120 


WEE Ba I 
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图 6 锡林郭勒 盟 地 区 积 雪 深 


ss， 即 亮 度 温 度 可 以 反映 
雪 内 部 的 信息 。 两 类 数据 在 积 雪 深度 估算 均 有 
rn, TERT 
射 率 数据 ( 表 1)。 本 文 首次 尝试 结合 两 类 遥感 数据 
用 于 积 雪 深度 的 估算 ,并 以 锡林郭勒 盟 e 
测 站 点 的 实测 数据 为 例 ,结果 表明 综合 利用 积 雪 表 
面 和 内 部 信息 明显 改进 了 积 雪 深度 估算 精度 。 
以 随机 森林 算法 为 代表 的 机 带 学 习 算 法 在 积 
雪 遥 感 领 域 的 应 用 越 来 越 广 , 这 些 学 习 算 法 能 够 处 
理 积 雪 深度 与 反射 率 和 亮度 温度 之 间 复 杂 的 非 线 
性 关系 。 机 器 学 习 算 法 (人 工 神 经 网 络 .支持 癌 量 
机 和 随机 森林 ) 均 可 以 取得 较为 满意 的 积 雪 深度 估 
算 精 度 ,尤其 随机 森林 算法 的 表现 最 好 。 尽 管 俩 最 
小 二 乘法 可 以 消除 变量 的 共 线 性 ,但 是 在 积 雪 深 度 
佑 算 领 域 的 表现 较 差 ,3 这 可 能 归 因 于 其 难以 适应 复 
杂 的 积 雪 场景 。 
本 文 提出 的 方法 综合 利用 了 光学 和 微波 遥 
的 优势 ,基于 机 央 学 习 算法 和 角 ee 
间 分 辩 率 的 积 雪 深度 空间 分 布 图 ,对 我 国 西北 地 区 


积 雪 深 度 遥 感 监 测 研究 具有 重要 意义 。 
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Estimation of snow depth based on reflectance and bright 
temperature in Xilin Gol League 


LI Hui-rong 
(Xilinhot National Climatological Observatory , Xilinhot 026000, Inner Mongolia, China) 


Abstract: Snow is an important water resource in arid and semi-arid regions of northwest China, which is also one 
of the important factors that affect the global climate change. Snow depth is the main parameter of the water 
resources and its estimation using a remote sensing technique has a vital effect on the assessment of resources and 
disease of northern China. The reflectance and brightness temperature data were the main data for snow remote 
sensing. Taking the snow depth data of 2012—2015 for Xilin Gol League, Inner Mongolia, China as an example, 
The MODIS/Terra satellite reflectance data and FY-3B brightness temperature data for the estimation of snow depth 
are combined for the first time in this study. The estimation accuracy for partial least squares and machine learning 
(artificial neural network, support vector machine and random forest) are then compared. When only reflectivity data 
is used, best accuracy is achieved by the random forest algorithm (RMSE=3.65 cm), followed by the neural network 
and support vector machine. The stochastic forest algorithm still performs best (RMSE=3.12 cm) when only the 
brightness and temperature data are used. However, no obvious difference is found with the artificial neural 
network. Using the partial least squares method exhibits the worst performance among the methods used. For the 
four methods, the accuracy of estimating snow depth using brightness temperature data is higher than that using 
reflectivity data, which may be attributed to the fact that reflectivity data is easy to saturate. Comprehensive 
research shows that the snow depth estimation accuracy based on the combination of reflectance and brightness 
temperature was better than that with single data. The random forest algorithm obtained the best performance, which 
is able to meet the needs of practical application. However, the accuracy of the three machine learning algorithms is 
not much different, with RMSE of 2.93-3.92 cm, which is far lower than that of the partial least squares method 
(RMSE=7.12 cm). Although partial least squares method can eliminate the collinearity of variables, its performance 
in snow depth estimation is poor, which may be attributed to its difficulty in adapting to complex snow scenes. 
Results of this paper have important significance for the study of water resources distribution and ecological 
environment assessment in northern China. 


Key words: snow depth; MODIS; MWRI; machine learning 


